Field-dependent Vortex Pinning Strength in a Periodic Array of Antidots 
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We explore the dynamic response of vortex lines in a Pb thin film with a periodic array of antidots 
by means of ac-susceptibility measurements. For low drive field amplitudes, within the Campbell 
regime, vortex motion is of intra-valley type and the penetration depth is related to the curvature 
of the pinning potential well, a. For dc-fields below the first matching field Hi, a reaches its highest 
value associated with a Mott Insulator-like phase where vortex lines are strongly localized at the 
pinning sites. For Hi < Hdc < H%, the response is mainly due to the interstitial vortices and a 
drops to smaller values as expected for this metallic-like regime. Strikingly, for H2 < Hdc < H3, 
we observe that a reduces further down. However, for Hz < Hdc < H4, a reentrance in the pinning 
strength is observed, due to a specific configuration of the flux line lattice which strongly restricts 
the mobility of vortices. We present a possible explanation for the measured a(Hd c ) dependence 
based on the different flux line lattice configurations. 



Aligned linear defects strongly influence the static and 
dynamic behavior of vortex lines in the mixed state. It 
has been extensively demonstrated that the introduction 
of columnar defects by heavy ions irradiation gives rise 
to several vortex phases and a considerable enhancement 
of the critical current. 1 ' 2 Theoretically, the expected vor- 
tex behavior should change significantly when the num- 
ber of vortices exceeds the number of columnar defects. 
However, the disorder in the randomly distributed pin- 
ning sites together with the pinning energy dispersion 
of the tracks lead to the appearance of matching ef- 
fects as washed out transitions or just barely defined 
crossovers. 3 In contrast to that, the introduction of na- 
noengineered periodic arrays of artificial pinning centers, 
where both topological and energetic disorder are absent, 
provides an ideal model system to explore the matching 
transitions and study the "zoology" of the vortex lattice 
configurations. 4 

In superconductors with periodic pinning arrays (PPA) 
the pinning properties are considerably improved at the 
applied field Hd c — Hi, where the density of vortices 
equals the density of pins. 4 This behavior is also present 
at higher commensurability fields H n — nHi, with n an 
integer number. However, the enhanced pinning progres- 
sively diminishes with increasing field since the vortex- 
vortex interaction becomes more relevant. An additional 
field dependence of the critical current arises due to re- 
arrangements of the vortex lattice configuration itself. 4 

Most of these vortex configurations have been di- 
rectly imaged by Lorentz and scanning Hall Probe 
microscopies, 5,6 and Bitter decoration 7 . However, none 
of these techniques allow one to extract the pinning po- 
tential depth. In order to do that, it is necessary to 
use indirect techniques like transport, dc-magnetization 
or ac-susceptibility. In transport and dc-magnetization 
measurements, the system is driven out of equilibrium, 
thus allowing the determination of critical currents and 



pinning strength of a moving flux line lattice (FLL), 
distinct of that observed by means of vortex imaging 
techniques. In contrast, low amplitude ac-susceptibility 
measurements permit one to probe the effective pinning 
potential without modifying the vortex distribution, thus 
making possible a direct correlation between FLL images 
and the observed dynamical response. 

Indeed, several theoretical and experimental works 8,9 
indicate that the ac response of a system with unidirec- 
tional defects exhibits, for small values of the alternating 
field drive h ac , a linear behavior characterized by an ac- 
susceptibility x independent of h ac and a very small dissi- 
pation. In this regime, known as the Campbell Regime, 10 
where vortices oscillate inside the pinning centers, it is 
possible to estimate quantitatively the strength of the 
pinning potential by analyzing the field dependence of 
the ac-susceptibility. 

In this work, we present measurements of the ac- 
susceptibility x = + °f a Pb thin film with a square 
periodic array of holes (antidots). We mainly concentrate 
on the linear regime where vortex displacement never 
surpasses the pinning range and the technique can be 
thought of as a passive non-destructive method to mea- 
sure the pinning potential. From the field dependence of 
the ac-penetration depth, A, we determine the curvature 
of the pinning potential well, a, for different matching 
fields. This allows us to estimate quantitatively, for the 
first time, the pinning potential for each specific vortex 
configuration. In particular, we show that the strength 
of the effective pinning potential exhibits a reentrant be- 
havior for i?3 < Hdc < Hi, where vortex mobility is 
strongly restricted. 

The studied sample was a 50-nm-thick Pb film of rect- 
angular shape (1.8 mm x 3 mm), with a square anti- 
dot lattice of period d — 1/im, which corresponds to 
the first matching field Hi = <&a/d 2 =20.7 Oe, with $ 
the superconducting flux quantum. The antidots have a 
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square shape with a size b — 0.25 /im. The maximum 
number n s of flux quanta that an antidot can hold is 
approximately given by n s « b/4£ (T), 11 where £(T) is 
the superconducting coherence length. For most of our 
samples we have estimated £(0) « 38nm and consider- 
ing the usual Ginzburg-Landau temperature dependence 
we obtain n s ~ 0.6 for T =6 K, the lowest temperature 
studied. Thus, we rule out the possibility of multiquanta 
occupation in a single hole. 

The film was electron beam evaporated in a molecular 
beam epitaxy apparatus onto liquid nitrogen cooled Si0 2 
substrates with a predefined resist-dot pattern and after 
the resist was removed in a lift-off procedure, the sample 
was covered with a protective Ge layer of 20 nm. 4 The 
sample has a T c =7.2 K with a transition width of 0.2 
K, for Hdc—0 and h ac —0.5 Oe. The ac-measurements 
were carried out in a commercial Quantum Design-PPMS 
device with drive field amplitudes h ac ranging from 88 
mOe to 10 Oe, and the frequency / from 10 Hz to 15 kHz. 
In all cases, the data were normalized by the same factor 
corresponding to a total step Ax' — 1, with Hd c =0. 

In order to identify the linear regime, we performed 
a series of measurements of x = x' + *x" as a func- 
tion of the ac excitation h ac for several temperatures at 
fixed / =3837 Hz and H dc =20.7 Oe. In Fig. 1 we show 
one of these measurements for T =6.4 K. The Campbell 
regime is characterized by a change in the magnetization 
dM = xhac proportional to h ac together with a scarce 
dissipation. We observe that x^ac departs from a linear 
response at h ac ~ 1 Oe. We also note that while in the 
linear regime the dissipation is negligible, an increase in 
the ac losses appears at the onset of nonlinearity. As we 
will show below, the range of ac-fields where x" ~ 0; is 
a reliable criteria to determine the range of the validity 
of the linear regime. 

An alternative way for determining the Campbell 
regime comes from the study of higher components of 
X- Indeed, nonlincarities should manifest themselves as 
higher harmonics Xm with m being an odd number. 12 We 
confirm this by measuring up to 9 components of x- We 
find that for the lowest h ac , Xm ~0 for m >1, while for 
h ac >1 Oe, x-m =/= for m=3,5,7 and 9. The acquisition 
of these components allows us to build the magnetization 
curve as a function of h ac for two extreme excitations h ac 
at T =7 K and Hd c =10 Oe, as is shown in the inset of 
Fig. 1. We clearly observe that for h ac =0.09 Oe the 
response is Meissner-like and almost free of dissipation 
(x" ~ 0). In contrast, for h ac =3 Oe the response can 
be satisfactorily described using the Bean Critical State 
model and the losses, given by the area of the loop, arc 
no longer small. According to the Bean model for a disk 
of radius R and thickness <5, the saturation magnetization 
M sat is related to the characteristic penetration field H p , 
through H p — 3M sat S/2R. 12 For the data showed in the 
inset of Fig. 1, 4irM sat - 1400 Oe and thus 2H p - 1.6 
Oe, in good agreement with the field needed to invert the 
critical state profile. 

Since the limits of the linear regime are determined 



mainly by the strength of the pinning sites, it is expected 
that the range in which this regime exists decreases with 
both, field and temperature, for a fixed h ac . This be- 
havior becomes evident in Fig. 2, where the out-of-phase 
component x" (upper panel) and the in-phase component 
x' (lower panel) are shown as a function of the scaled field 
Hdc/ Hi for two values of h ac at T =7 K. 

For small amplitudes (h ac =0.09 Oe), and for < 
Hdc < H± a high screening and very small losses are 
observed. However, for Hd c > Hi, the screening exhibits 
an overall decrease together with a sequence of minima 
at Hdc = H n caused by local depressions in the losses. 
This behavior can be understood within a scenario based 
on the picture of vortices oscillating inside their pinning 
sites, without large excursions. Indeed, at the commen- 
surability fields, where the effective pinning is maximum, 
a higher screening and a reduction of dissipation due to 
the more confined movement of vortices, is expected. 

On the other hand, at h ac =1 Oe, deep within the crit- 
ical state, the response is remarkably different. Although 
the local dips at H dc = H n in the screening curve are still 
present, now these minima are accompanied by a corre- 
sponding peak in the losses. This is a fingerprint of the 
critical state regime. 13 At this high amplitude value, vor- 
tices are no longer localized in small pinning valleys, but 
instead they can travel larger distances (inter-valley mo- 
tion). As we pointed out above, in this regime, losses are 
given by the area of the M/i ac -hysteresis loop (see inset of 
Fig. 1) whose width is determined by the critical current 
of the system. At H dc = H n , the critical current reaches 
a maximum and therefore also the losses are maximized, 
as is indeed observed in our experiments. Based on the 
previous analysis we will adopt as a conservative and safe 
criterion to ensure a linear regime, the limit of the used 
Hdc range, to the zone where x" <0.01. 

In the upper panel of Fig. 2 we have also included 
x"{H dc ) curves for T =6.8 K and 6.5 K, at h ac =0.09 Oe. 
As expected, we find that with decreasing T, the range 
where the system is free of dissipation, becomes broader. 
In particular, at T =6.5 K, the losses are negligible for 
the whole field range studied. 

From now on, we will concentrate on the field depen- 
dence of x m the linear regime for h ac =0.09 Oe, at 
two different temperatures, T =6.5 K and 6.8 K. In this 
regime, the ac magnetic response is determined by the 
complex penetration depth X(B ~ Hd Cl T, /) = Xr + iA/. 
For ideal systems without dissipation, A is real and 
x" = 0. However, dissipation leads to a / / and 
a ratio e = Xi/Xr, although small, is finite. According 
to ref. 14 , x' is solely determined by A and the sample 
geometry through x ~ -1 + X^°[ c «/( A n + ¥>)]> where 
ip = C/X 2 , C is a constant depending on the sample 
shape, and A„ and c n are real numbers. We can invert 
this equation in order to compute A from our x data, as 
described in ref. 15 . 

For our experimental conditions e <C 1 and A 2 w 
X 2 R (B,T) = X 2 L {T) + X 2 C {B,T), where X L {T) and 
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\c(B,T) = ^/B$>o/4:ira are the London and Campbell 
penetration depths, respectively. The Labusch constant, 
a, represents the curvature of the pinning well assumed 
to be parabolic for small displacements. 

Fig. 3 shows \ 2 R as a function of the reduced field 
Hd c /Hi for T =6.5 K. The most evident feature of this 
figure is that for well defined field ranges (0 < H n -\ < 
Hdc < H n < H 5 ), \ 2 R follows a linear field dependence, 
\ 2 R oc B/a n , in agreement with the Campbell regime 
scenario, with a Labusch parameter a n clearly depend- 
ing on the field range. This is one of the main results 
of this work. The same behavior is observed at T =6.8 
K (see upper inset of Fig. 3) although, as we previously 
noticed, the linear regime is, in this case, constrained to 
Hdc < H 3 . 

Since at H dc =0, \ 2 R = A|(T) = A|(0)/(1 - t), where 
t = T/T c , using the data of Fig. 3 we can easily estimate 
the London penetration depth at T =0. The absolute 
value Al(O) «52 nm so obtained, turns out to be rela- 
tively close to the accepted value Al(0) w 47 nm. 16 

In all cases, A < and therefore the penetration of 

the compressional vortex waves involves just a few units 
cells of the antidots array from the sample border. As 
soon as a new vortex configuration develops right after 
a commensurability field H n , with a consequent change 
in the average pinning strength, it will be simultaneously 
detected in the ac susceptibility response. Since our field 
step is 1 Oe, the density of the new developing FLL con- 
figuration is high enough to record a sharp change at H n . 

How actually the new vortex configuration penetrate in 
a sample with PPA has been a matter of controversy dur- 
ing the last years. 17,18 Cooley and Grishin 17 have shown 
that in a one-dimensional system with PPA, a terraced 
flux profiles should be established. In our system (thin 
films), due to the small gradient in the flux profile, a sin- 
gle terrace critical state is expected. According to this 
model, increasing the field beyond H n , leads to the mo- 
tion of the new vortex configuration towards the center of 
the sample. Since we are sensing the magnetic response 
at the sample border, the same A 2 should be measured as 
Hdc increases and therefore no field dependence should 
be observed. This is in sharp contrast with our observa- 
tions. Instead of that, we argue that for each applied field 
a mixture of the patterns corresponding to the nearby 
matching fields distributes in the whole sample. 18 Within 
this scenario, a so obtained results from an average of the 
different pinning potentials associated with each existing 
vortex configurations. 

In the lower inset of Fig. 3 we show a zoom of the same 
data as in the main panel, for H dc < Hi. We observe 
that for Hdc < Hi, where vortices are strongly pinned 
by the antidots, a is larger than for Hd c > Hi where 
caged interstitial vortices, much weaker pinned and with 
higher mobility, dominate the ac response. The Labusch 
parameter obtained from these slopes is summarized in 
Fig. 4 for all fields studied. The sketches of different 
stable vortex configuration according to previous vortex 
imaging experiments 5 and simulations, 19 are also shown. 



Since for Hd c > Hi the entity which respond to the 
external ac excitation is no longer a well identified single 
vortex, but a complex movement involving all interstitial 
vortices in a unit cell, the interpretation of the ac sus- 
ceptibility for these higher fields becomes more puzzling. 
For instance, at Hd c = Hi, vortices start to form inter- 
stitial dimers which can adopt two different orientations 
separated by a small energetic barrier. When the drive 
field is applied, a librational motion of the vortex pairs 
is excited. Since the observed Labusch parameter is even 
smaller than for Hd c < Hi we argue that this motion 
permits (in average) larger excursions of the vortex lines 
than the movement of single interstitial vortices at the 
same h ac . 

For fields higher than H$ , the unit cell motif consisting 
of a rhombic configuration, loses the rotational degree 
of freedom and consequently, a reentrance of a(Hd c ) is 
observed. On top of that, because of the system is now 
changing from a phase with higher mobility to another 
phase with lower mobility, the density of the new unit 
cell motif needed for significantly change the measured 
a, should be larger. As a consequence, the transition 
does not occur exactly at H% but a field slightly higher, 
as we indeed observe. 

For fields H dc > H4, vortices in a unit cell distribute 
forming a square arrangement where the rotational de- 
gree of freedom is recovered, thus giving rise to an easy 
librational movement and a lower average pinning poten- 
tial. At fields higher than the fifth matching condition 
nonlinearities develop and no reliable Labusch constants 
can be determined. It is important to stress that the 
linear dependence of \ R with field for Hd c > Hi is a 
highly non trivial result. Clearly, the complex mecha- 
nisms involved during the ac excitation deserves further 
experimental and theoretical investigation. 

In general, the main trends shown in Fig. 4 are similar 
to those reported in simulation studies made recently by 
Reichhardt et al. 20 In that work, the authors determine 
the dependence of the critical depinning force, as a func- 
tion of H dc / Hi for a system with a square pinning array. 
In agreement with our experimental observations they 
find a series of roughly constant plateaus of decreasing 
value. However, unlike the Labusch constant, the de- 
pinning force density and the critical current depend on 
the depth of the pinning potential (rather than its cur- 
vature) as well as the FLL configuration, and therefore 
this analogy should not be pushed too far. Essentially 
transport and ac-measurements in the linear regime are 
not directly comparable. 

Finally, it is worth noting that in similar studies per- 
formed in high temperature superconductors with colum- 
nar defects introduced by irradiation, 9 the change in the 
slope at the first matching condition turned out to be a 
blur crossover with an onset around Hi/ 2, presumably 
due to the topological and energetic disorder of these 
defects. 3 In those studies a field dependent Labusch pa- 
rameter was found for Hd c > H\/2, as a result of both, 
a disordered state where a single unit cell motif can not 
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be identified, and the lattice interactions when the sys- 
tem approaches to the collective regime. In our system 
with PPA where vortices of same species have the same 
surrounding, the regimes are connected by abrupt tran- 
sitions and a remains constant for the low fields studied. 

In summary, we have used the ac-susceptibility tech- 
nique as a non-disturbing tool for determining the influ- 
ence of the vortex lattice arrangement in the averaged 
pinning strength. We have estimated quantitatively the 
value of these effective pinning potentials for several sta- 
ble vortex configurations. This powerful method has the 
main advantage of an easy extension to different systems 
with periodic pinning array and the possibility of a direct 
comparison of the pinning energies found from different 
experiments. 
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FIG. 1. Drive dependence of dM = xhac at the first match- 
ing field Hi, T =6.4 K and / =3837 Hz. The arrows indicate 
the onset of the nonlinear regime. The inset shows the nor- 
malized magnetization derived from several harmonics Xn for 
two extreme amplitudes h ac =0.09 Oe (filled symbols) and 
h ac =3 Oe (open symbols), at T =7 K. 

FIG. 2. Field dependence of the ac-susceptibility 
X = x' + *x" components for h ac =0.09 Oe (filled circles) 
and hac =1 Oe (open circles) at T =7 K. In the upper panel 
is also shown x" f° r T =6.5 and 6.8 K, at h ac =0.09 Oe. 

FIG. 3. Field dependence of \ R at T =6.5 K (main panel) 
and T =6.8 K (upper left inset), for /i ac =0.09 Oe. The lower 
right inset shows a low field zoom of the main panel data. 
The horizontal and vertical units in the insets are the same 
used in the main panel. The lines are guide to the eye. 

FIG. 4. Field dependence of the Labusch parameter a 
in units of its maximum value a max at T =6.5 K. The 
sketches represent the different stable vortex configurations 
where empty antidots are drawn as open circles, occupied 
antidots as big filled circles and interstitial vortices as small 
filled circles. 
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